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Abstract. The idealized model of the geomagnetic field-solar wind interaction yields a 


singular (neutral) point on the magnetopause at which the magnetic field vanishes. If we 
expand the fields in a power series around this point, including quadratic terms, we can derive 
an approximate equation of the magnetopause in a small neighborhood of the neutral point 
w’hich is consistent with the idealized boundary conditions to fourth order. We then consider 
the additional pressure due to multiple reflections of particles in this neighborhood and show' 
that less than a 4 per cent correction to the single-reflection pressure condition is necessary. 


Introduction. In the idealized model of the 
geomagnetic field-solar wind interaction [Beard, 
1960] we assume the field to be excluded from 
the plasma and contained in a cavity called the 
magnetosphere. The boundary of the cavity, or 
magnetopause, is unknown, but we stipulate the 
dynamic condition that the magnetic pressure 
just inside the magnetopause be exactly balanced 
by the kinetic pressure of solar wind particles 
elastically reflected from the surface (thermal 
effects are neglected). Figure 1 is a schematic 
drawing of the magnetosphere with the dipole 
located at the origin of coordinates and per- 
pendicular to the direction of the plasma stream. 

Mathematically, this situation is described by 
the field 

H = —grad 0 V 2 fl = 0 (1) 

inside the magnetosphere, and the boundary 
conditions 

H-grad F = 0 (confinement) , , 
H-H = (3 2 cos 2 x (pressure) 


gradient is in the x' direction. Although the 
problem is unsolved as yet, there are some good 
approximations to the general shape [Mead and 
Beard, 1963; Spreiter and Briggs, 1962; Midgley 
and Danis, 1963]. However, these approximations 
generally fail in the region of the neutral point, 
a region of interest since, if it does exist, it is 
likely to be unstable, and a possible point of 
entry for high-energy particles into the mag- 
netosphere. 

The neutral point considered here is an X-type 
neutral point [Dungey, 1958, pp. 39-41, 51-52, 
98-102]; however, it lies on a bounding surface, 
w'hich must ran parallel to one of the limiting 
field lines (Figure 2). We must assume that 
surface currents can be made to account for the 
disappearance of the field outside the surface 
(magnetopause) . 

Series representation. We can represent the 
field H = /3h by means of its scalar potential S2. 
We first transform to the ( x , y, z ) coordinate 
system with origin at N such that 

R(dx, dry, dz) = (dx r , dy' , dz') (3) 


which hold on the surface, F{x', y', z') = 0, 
where x is the angle of incidence of the incoming 
plasma (Figure 1), and /3 2 = 8Tr(2nmV 2 ). V, n, 
and m are the plasma drift velocity, ion number 
density, and ion mass, respectively. 

A singular point on this unknown surface is N, 
the neutral point, where the magnetic field lines 
‘split,’ the field vanishes, x = w/2, and the 
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In a region small compared to the apex radius R, 
w r e have x* + y 1 + z 2 <<C 1. The scalar potential 
is expanded in the form 

— Q./B = <tx 2 + by 2 + cz 2 + dx + ey s + jz 

+ mxy 2 + nxz + rzx + szy 2 (4) 

We must include cubic terms in the potential, 
since w r e expect the limiting field lines to be 
curvilinear in the xz plane. There are no linear 
terms, because the field must vanish at N ; terms 
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and the requirement V • H = 0 yields 
a b d~ c = 0 3d m n — 0 



y 

Fig. 1. Schematic drawing of the magnetosphere. 


in y to the first power are deleted, because 
symmetry requires H, = 0 in the xz plane. 
The dependence on xz has also been omitted, 
since for y = 0 = z there should only be an 
x component of H on the x axis (Figure 2). 

The gradient of (4) yields the normalized field: 

h z = 2ax + 3 dx 2 + my 3 + nz + 2 rxz 

A y = 2by + 3 ey 2 + 2 mxy + 2 szy (5) 

h, = 2 cz + 3/z 2 + 2 nxz + rx + sy 2 


( 6 ) 

3/ + r + s = 0 e = 0 


Furthermore, as Figure 2 shows, we expect 
a < 0, c > 0. 


Noon meridian contour. In the 

XZ 

plane the 

boundary conditions (equation 2) 

on 

the 

noon 

meridian contour become ( dx/dz = 




K/h. = -FJF Z = 

x' 


(7) 

h z + A, 2 = x' 2 /(l + 

z /2 ) 


(8) 

If we substitute for x' in (7) and (8), ' 

then 


hj + h 2 (2h, 2 — 1) + h 

\ 

0 

(9) 


and, since both components must vanish at N, 
h? - (1 - 2 h, 2 - (1 - 4/i, 2 ) 1/2 )/2— > h* 


N (10) 

From Figure 2 we see that h x is opposed in sign 
to the coordinate z on the noon meridian contour. 
Thus, sufficiently close to N, 

h x = ±/i, 2 for 2^0 (11) 

If we substitute into (7) and differentiate, 

x" — ±dh,/<k z 5 0 (12) 


X 


Fig. 2. 


Fields near the neutral point in the noon meridian plane. 
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which leaves two alternatives: either x" is 
discontinuous at N or x" = 0 at N. The latter 
implies c = 0, which is inconsistent with the 
geometry of the limiting field line in the xz 
plane. Thus the curve, {x = x(z), y = 0}, must 
be represented as two separate power series for 
z ^ 0. Since (11) must be satisfied on the noon 
meridian contour, it gives an implicit representa- 
tion of that contour near N. This must agree 
with x = / dzhjh , of (7) and with Figure 2. 
The leading term of (11) is 

2 ax + nz = ±4c 2 z 2 (13) 

If n = 0, equation 13 implies that a is of positive 
sign, which contradicts Figure 2. Therefore we 
must assume that n has two different values, 
according to the sign of z. This is permissible as 
long as h remains continuous at z = 0, and 
V-h = 0 everywhere. If we compare (13) with 
the leading term of (7), we find 

x' — ±h t <0 x — ±cz (14) 

near N. Thus, (13) becomes 

n — db 2c(2c — a) for z ^ 0 (15) 

Surface representation. Since the gradient is 
in the x direction at N, we can represent the 
surface by two second-order expansions in (y, z) : 

F(x, y, z) 

= x - AV - B*y 2 - C*yz = 0 (16) 

for z ^ 0. By symmetry, B + — B~ = B; (14) 
requires A + = — A~ = c. If we form the dot 
product H-VF = 0 and substitute for x from 
(16) we find the confinement condition satisfied 
to second order in z 8 , y i , and yz when 

m = 2B(2b — a ) (17) 

C*(l - b/a - c/a) = 0 (18) 

The divergence conditions (equations 6) require 
that d also be discontinuous, since 3d -j- 
m + n = 0, although m is continuous. There- 
fore, from (17) and (15), 

3d = ±2c (a - 2c) + 2 B(a - 26) z ^ 0 (19) 

Equation 18 is incompatible with a + 6 + c = 0 
unless C* = 0. Although (19) forces h x to be 
discontinuous at z = 0 (o^ 2c), the 3 dx* term 
is negligibly small compared to the ax term in (5). 
Third-order terms in H-VF = 0 are elimi- 


nated by simply setting / = 0 = s. Then, by 
the divergence condition, r = 0 also. Thus in 
the power series representation, the confinement 
condition of (2) is satisfied up to fourth-order 
errors, on the postulated surface of (16). If we 
substitute the expansions for h and F(x, y, z) 
into the pressure condition of (2), we find 

hj 3 + h 3 + h, 2 = cos 2 x 

= {dF/dzf{\ — 2nd-order terms) (20) 

and expanding h yields 

4c 2 z 2 + 4 }/ y 2 -J- 4th- order terms 

= 4cV(l - 2nd-order terms) (21) 

By setting 6 = 0 we can reduce the error in the 
second boundary condition to fourth order. Our 
final expansion for the field in the neighborhood 
of N is given by : 

— 0//8 = cz — cx — (m ± 6c 2 ) x / 3 

+ mxy rfc 6c 2 xz 2 (22) 

and the surface of the magnetosphere near N is 
given by 

F(x, y,z) = x — (±)cz 2 — my 2 /2c (23) 

for z ^ 0. 

Multiple reflections at N. Midgley and Davis 
[1963] have observed that the effect of multiple 
reflections of particles near N might seriously 
alter the pressure condition (equation 2) near 
the neutral point. In the second-order approxi- 
mation, we can show that the shape of the 
surface is consistent with the pressure condition 
to within 4 per cent. This effect would be most 
pronounced in the xz plane; we shall calculate 
the added pressure bp at a point (x, z) on the 
noon meridian contour due to multiple reflections. 

Figure 3 represents a particle incident on the 
point (lo, z«) on the contour x = cz 8 , z < 0, 
y = 0, with an angle of incidence Xo- This parti- 
cle, on first reflection, strikes the magnetopause 
again at (x, z) where 

(x — x 0 )/(x — 2o) = c(z + z 0 ) = tan 2x 0 (24) 

and 

cot xo = — x'(z 0 ) = — 2cz 0 (25) 

therefore, 

z = 4cz 0 /(l — 4c 2 z 0 2 ) — z 0 = 3z 0 (26) 
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as indicated in Figure 3. We note that the angle 
of incidence is ir — 2xo + x at (z, z), and the 
change in the normal component of momentum 
will be correspondingly smaller. If we pursue the 
particle to the third point of reflection we find 

c(z + 3z 0 ) = tan 2(x — Xo) z — 5z 0 (27) 

However, the effect of the particle, by this time, 
is negligible. 

We now consider the additional pressure bp 
at {x, z) due to particles reflected from (x 0 , z 0 ) 
and compare it with the pressure p = 2 nmV* 
cos* x of the incident plasma stream (ignoring 
the y coordinate). If nV dx 0 particles/sec are 
incident on an area of ds 0 = (dxo 1 + <fz 0 ! ) U2 at 
(x„, z 0 ), they are reflected onto an area ds = 
(<2x 2 + dz 2 ) 1,J at (x, z), where their momentum 
changes by 2 mV cos (ir — 2xo + x) per particle. 
Therefore, since dx = 9 dx 0 , 

8p = —2 mnV 2 cos (x — 2xo)(dx/ds)/9 (28) 

If we substitute for x> Xo from (25) and note 
that dx/ds = cos x> then we find 


8p/p = 1/ 27 « 1 (29) 

We can conclude that the effect of multiple 
reflections near the neutral point is negligible 
and that the pressure condition is valid over the 
entire magnetosphere in this idealized repre- 
sentation. 
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